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Abstract 



Over the last few years, the fruitful data provided by the Large Area Telescope aboard the Fermi Gamma- 
ray Space Telescope has revolutionized our understanding of high-energy processes in globular clusters, 
particularly those involving compact objects like millisecond pulsars (MSPs). Gamma-ray emission be- 
tween 100 MeV to 10 GeV has been detected from more than a dozen globular clusters in our Galaxy, 
most notably 47 Tucanae and Terzan 5. Based on a sample of known gamma-ray globular clusters, 
empirical relations between the gamma-ray luminosity and properties of globular clusters such as stellar 
encounter rate, metallicity, as well as optical and infrared photon energy density in the cluster, have 
been derived. The gamma-ray spectra are generally described by a power law with a cut-off at a few 
GeV. Together with the detection of pulsed 7-rays from a millisecond pulsar in a globular cluster, such 
spectral signature gives support that 7-rays from globular clusters are collective curvature emission from 
magnetospheres of MSPs within the cluster. Alternative models in which the inverse-Compton emission 
> ■ of relativistic electrons accelerated close to MSPs or the pulsar wind nebula shocks have also been sug- 
gested. Observations at >10 GeV by Fermi/LAT and atmospheric Cherenkov telescopes like H.E.S.S.-II, 
MAGIC-II, VERITAS, and CTA will help to settle some questions unanswered by current data. We also 
discuss TeV observations of globular clusters, as well as observational prospects of gravitational waves 
1 from double neutron stars in globular clusters. 
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Globular clusters are the oldest gr avitationally- bounded stellar systems in the Galaxy. Nearly 160 



globular clusters are known nowadays (|Harrist 1 19961 2010 edition). They fill a spherical halo around the 



Galaxy, many of them are located within the Galactic bulge. Due to the high concentration of stars 
within the globular clusters, they host a large number of compact objects including neutron stars and 
white dwarfs; many of them are found in binary systems, forming for example low-mass X-ray binaries 
(LMXBs) and Cataclysmic variables. 

Since 1970s, it has been known that the formation rate per unit mass of LMXBs (Alpar et al. 1982) is 
orders of magnitude greater in globular clusters than in the rest of the Galaxy (Katz 1975; Clark 1975). 
As millisecond pulsars (MSPs) are generally believed to be descendants of LMXBs, it becomes natural 
that ^80% of the known MSPs are detected in globular clusters (cf. Manchester et al. 2005). Theoretical 
arguments have long asserted that the formation of LMXBs (and therefore their decedents MSPs) is 
made efficient through frequent stellar encounters. Using the X-ray populations in various globular 
clusters unveiled by the Chandra X-Ray Observatory, Pooley et al. (2003) and Gendre et al. (2003) 
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found a positive correlation between the number of LMXBs in globular clusters and the stellar encounter 
rate, r c , putting the dynamical formation scenario of LMXBs in globular clusters on an observational 
ground. 

Using the cumulative luminosity distribution functi ons of radio millis econd pulsars (MSPs) in globular 



Hui et al 



(2010) found that the number of MSPs 



clusters as a probe of the resided MSPs in the clusters, 
in a globular cluster is correlated with its stellar encounter rate, as well as its metallicity. This finding 
provides the first observational evidence of the dynamical origin of MSPs. This is easy to understand 
since MSPs are descendants of LMXBs. 

In this paper we review main results revealed by high-energy 7-ray observations of Galactic globular 
clusters, which were mainly contributed by data taken using the Large Area Telescope (LAT) on board 
the Fermi Gamma-ray Space Telescope since its launch in June 2008. Towards the end of the paper we 
also discuss the recent TeV observations as well as observational prospects relate d to gravit a tiona l wave 
experiments of globular clusters. The readers are referred to an earlier review bv iBednareld ([20111 ) from 
a more theoretical aspect. 



2. GeV 7-ray observations 



2.1. Pre-Fermi era 



Observations by EGRET in the 1990s did not result in any positive detection from a c ouple of globular 



clust e rs, including 47 Tucanae that hosts the largest number of known MSPs at that time (jMichelson et al 



1994; 



Manandhar et al 



19961) . Later LAT observations found that the EGRET upper limit for 47 Tu- 



canae is only a factor of two higher than the measured flux. 



2.2. First discoveries: Jf.1 Tucanae and Terzan 5 

The launch of the Fermi Gamma-ray Space Telescope has enabled the detection of 7-ray emission 
from a numb er of globular cluste rs. The first discovery is made on 47 Tucanae by the Fermi/LAT 



collaboration (jAbdo et all l2009bl see Fig. []}. By th e time this discover y was announced, MSPs as a 



class had already been established as 7-ray emitters (jAbdo et all 12009a) , and they are believed to be 



the only stable 7-ray sources in globular clusters. T he observed 7-rays f rom the direction of 47 Tucanae 



are naturally attributed to the MSPs in the cluster ([Abdo et all l2009bf ). This assertion is strengthened 



2009c 


). From the 


s" 1 ( 


Abdo et al.. 



averag e 



they woul d not be det e cted at distances of several kpc (the very luminous PSR J1823— 3021A is an 

see Sect. 



exception; 



Freire et al 



2011 



It is therefore generally believed that the 7-rays from 
globular clusters do not come from a single MSP, but are collective emission originated from the entire 
MSP population in a cluster. Assuming that the average 7-ray efficiencies of MSPs in 47 Tucanae are 
similar to nearby MSPs, the number of MSPs n eeded to give rise to the observed flux from 47 Tucanae 
(-2.5 xlO" 11 erg cm" 2 s _1 (jAbdo et all boiObh l is about 50, a factor of -2 higher than 23 MSPs in 
47 Tucanae detected in radio survey. Such discrepancy can be understood as incomplete radio surveys of 
MSPs, e.g., due to viewing angle effect. 



The second globular cluster that is found to emit 7-rays is Terzan 5 ([Kong et all l2010l see Fig. [T]). 
Terzan 5 hosts the largest number of detected MSPs by now (i.e., 33). Its spectrum also shows cutoff 
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Figure 1 : Photon count maps of 47 Tucanae (up per left ; Ubdo et all l2009bh and Terzan 5 (upper right; iKong et all bold) . 
Their 7-ray spectra are also shown lAbdo et allboiObh . 
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100 MeV is given by -7.1 xl0~ n erg cm ~ 2 s~ 1 (lAbdo et al 



47 Tucanae and Terzan 5 to be 4.0 kpc ([McLaughlin et al 



which is at —3 GeV (jKong et all l2010t lAbdo et al. . 2010b , s ee Fig . HR). The observed energy flux above 



2010b) l^j. Taking th e distance estimates o f 
b006h and 5.5 kpc (|Ortolani et all l2007h . 



respectively, the 7-ray luminosity of Terzan 5 is roughly five times that of 47 Tucanae, indicating that 
either the actual number of MSPs is higher in Terzan 5 than 47 Tucanae, or the average 7-ray efficiency 
of the underlying MSPs in Terzan 5 is higher than 47 Tucanae, or a combination of both. There are some 
indications from c umulative radio luminosity distributions that Terzan 5 may indeed contain more MSPs 
than 47 Tucanae ( Hui et all I2010I) . If the number of MSPs resided in Terzan 5 is smaller than —5 times 



that of 47 Tucanae, environmental factors such as the bac kground soft photons (whos e nu mber is much 



higher in Terzan 5 than in 47 Tucanae) , as suggested by iBednarek fc Sitarekl (|2007f ) and I Cheng et al 
(2010) , may explain t he higher 7-ray luminosity of Terzan 5. 



Kong et al 



(2010) report a marginal detection (3.7a) for 7-ray emission between 10 GeV and 20 GcV 
for Terzan 5 and no sign for emission in the same energy band for 47 Tucanae. However, this may not 
necessarily indicate a substantial difference between 47 Tucanae and Terzan 5 in their >10 GeV spectra. 
The energy flux of Terzan 5 is higher than 47 Tucanae by a factor of —3, which might explain the apparent 
discrepancies in the 10-20 GeV band for these two globular clusters. 

2.3. Globular clusters as a population of "/-ray sources 



After these two initial discoveries. 



Abdo et al 



Tarn et al 



(|2011r ) have reported detailed 



(l2010bh and 

analyzes of a total of 16 globular clusters. These include 12 cases for which the 7-ray detection is sig- 
nificant: Terzan 5, Liller 1, 47 Tucanae, NGC 6440, M 62, M 28, NGC 6388, w Centaurus, NGC 6624, 
NGC 6652, NGC 6441, and NGC 6752 (in descending order of reported >100MeV energy flux in the 
best-fit spectral model), and four cases for which the 7-ray detection is not significant enough yet: 
NGC 61 39, NGC 6541, M 15, and M 80. There are two more detections in the second Fcrmi/LAT 
catalog ([Nolan et alll2012l ): IC 1257 and 2MS-GC01. The 7-ray properties for these 18 globular clus- 
ters are shown in Table [1] No variability has been found for any reported detection, indicating that the 
7-ray emission is stable. 

In the original analyzes, for almost all of the cases (except Liller 1 and NGC 6624) where the 7- 
ray emission is significant enoug h, the nominal clust er position is within i?gs from the best-fit 7-ray cen- 
troid positions. For NGC 6624, 



Freire et al 



([20 111 ) has refined the 7-ray position, which is consistent 
with the core of the cluster. In the case of Liller 1, no strong radio and X-ray counterpart is found 
near the 7-ray emission region and therefore the MSP population in the respective clusters are the only 
known potential source of the observed 7-ray emission. In the context of magnetospheric origin of 7-rays, 
the offset requires that most MSPs in these clusters are located outside the core, a possibility which is 
unlik ely. In the inverse- Compton model, relat i vistic electrons diffuse gradually outward from the cluster 



core ( Bednarek fc Sitarek 



2007 



Cheng et al 



20101 ) . If the diffusion time scale of the accelerated parti- 



cles is much shorter than the cooling time scale, the emission from these particles may be offset from the 
clu ster cores if t hey dif fuse unevenly into different directions. 



Abdo et al 



(|2010bf) report that the 7-ray emission region is consistent with a point source for almost 
all globular clusters that they analyzed, except possibly M 28. The upper limit of the angular size of the 
emission ranges from 4.8' (for 47 Tucanae) to 15.6' (for M 28). 



1 The energy flux over 0.1-10 GeV derived in iKong et alj J2OI0I) 



^1.2xl0 -11 erg cm _2 s 



which is a factor of ~2 



higher. This may be due to the differences in the choices of region-of-interest (ROI) and energy range in the two publications. 
Here we quote the results in lAbdo et al.l fcoiObh for both 47 Tiicanae and Terzan 5 to better compare their 7-ray properties. 



Table 1: The 18 reporte d 7-ray emitting glo b ular clusters and candidates (whose 7 -ray detection significance <7.0<r) until 
end of 2011. References: Ubdo et all i hoiObh ; IrTolan et ail j2012h ; iTam et al.l fedllh 
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The spectra of the 7-ray emission from globular clusters generally show signs of spectral cut-off at 
energies 1-4 GeV, with different significant levels (see Table[lJ. On the other hand, possible enhancements 
at high energies (up to ~40 GeV) are reported for Liller 1 and NGC 6624, although in the latter case 
there may be contamination by a nearby flaring radio source, detected at 4-ct level in the first year of 
data. Future observations at energies >10 GeV are crucial to probe or rule out a possible enhancement 
at high energies for globular clusters, thus shedding light on the origin of 7-rays from globular clusters. 



2.4- Searching for j-ray pulses from globular clusters 

Attempts have been made to detect pulsed 7-rays from individual MSP in globular clusters. In the 
first three years after the launch of Fermi, no 7-ray pulsation was found using available ephemerides for 



MSPs in various globular clusters (jAbdo et al 



2009b 



2010b ) 



The first successful effor t to detect 7-ray pulsation from an individual MSP in a globular cluster was 
(|2011l ). They found a significant, 7-er detection of 7-ray pulsations above 100 MeV 



Freire et al 



made by 

from pulsar J1823-3021A in the globular cluster NGC 6624 (Fig. 0). The high measured flux in the 
pulsed component significantly constrain the total number of MSPs in NGC 6624 to be smaller than 32 
and provide a strong support to the pulsar magnetospheric emission model of 7-rays at least from this 
globular cluster. 



M 28 hosts 12 MSPs, the third largest among known GCs (jRansoml . I2008T) . after Terzan 5 and 
47 Tucanae. The cluster contains the MS P, M28A/PSR J18 24-2452A, that is very energetic and is the 
first MSP discovered in a globular cluster (jLvne et all Il987l ). The AGILE collaboration claimed a >4-er 

5 



hint of pulsation at >100 MeV during only a 6-day period in 2007 (IPellizzoni et all 120091) . However, this 
tentative evidence has not been confir med by the Fermi/LA T collaboration, who found no pulsation from 
PSR J1824-2452A since their launch (|Abdo et all boiObh . 



3. Which properties of globular clusters determine the 7-ray luminosity? 



3.1. Properties of globular clusters related to the j-ray emission 

The two-body stellar encounter rate has lo ng been believed to relate to the bin ary formation rate and 



hence the number of MSPs in a GC (see, e.g. iPoolev et al 



2003 



Huietal 



201011 . The first at t empt t o 



search for correlation of 7-ray luminosity with the encounter rate was preformed by 



Abdoetal 



(l2010bh . 



Using the first 8 known 7-ray globular cl usters and the express ion T 



PC) 



where po is the central 



luminosity density and r c the core radius (jVerbunt fc Hutl 119871 ). they found that the 7-ray luminosity 
of a globular cluster is correlated to the encounter rate with a linear correlation coefficient 0.7. Assuming 
that each MSP emits similar amount of 7-rays, the number of MSPs in t hese 7-ray emitting globular 
clu sters is said to b e consistent with the number given by other estimators (jAbdo et alll2010bl ). 

(1201 II) used the expressio n T c oc p\r\j<j where a is the velocity dispersion at the cluster 



Hui et al 



center calibrated by 



Gnedin et al 



(|2002D . Using the 7-ray data of 15 globular clusters, they found that the 
7-ray luminosity of globular clusters is correlated to the encounter rate at 99% level, as oc T c 

Besides the stellar encounter rate, there are oth er parameters that are thought to relate to the number 
of MSPs in a globular cluster. 



Hui et al 



(2010) identified the metallicity of a globular cluster to be 
another important indicator to the number of MSPs. A binary system in a globular cluster with a 
higher metallicity is more efficient in orbital sh rinkage by mag netic braking. This gives rise to a higher 
likelihood of a successfu l Roch e-lobe overflow (jlvanoval . 120061 ) and in turn leads to a higher formation 
rate of MSPs. 



Hui et al 



(|2010l ) found a positive correlat ion between met allicity and the MSP population 



Using the available 7-ray data from 15 globular clusters. 



Hui et al 



([201 II ) found that the 7-ray luminosity 
of globular clusters is correlated to the metallicity, [Fe/H], at 99.9% level, as L 1 oc r0.59±o.i5^ 

On the other hand, primordially-formed binaries are not related to stellar encounter rate. If they 
form the majority of binaries, one would expect the binary population to be correlated to the cluster 
mass, Mqc- Assuming a constant mass-to- light ratio, Mqc can be estimated from the absolute visual 
magnitude My according to Mqc = 10~ 4Mv . This mass estimate has the advantage of not having any 
correlation with the encounter rate, thus it serves as an independent parameter in the correlation study. 
However, no correlation of My with the 7-ray luminosity was found (jHui et all l201lh . 

Apart from the above three cluster parameters, the energy densities of optical, u opt , and infrared 
photons, ujr, at the cluster location are also important in t he predictions of the 7-ray properties as they 
serve as the seed photons in the inverse Compton models (jCheng et all l2010f) . Using the GALPROP 
code (jStrong fc MoskalenkoL I1998T) to estimate the soft photon energy densities, the 7-ray lu minosity is 



found to correlate with them at >96% level, as L- 



oc Uop 7 t 8± °' 27 as well as L 7 



,1.29±0.44 
X IR 



( Hui et al 



20111 ). While the original study was based on the Harris ca talog (2003 v ersion), we have re-examined the 
above correlations with in the revised 2010 Harris catalog (jHarrisl . 119961 ) and found that the correlations 
remained robust. 



3.2. Fundamental planes of ^ -ray emission from globular clusters 

We have shown that the total energy output in 7-rays indeed scales with the two factors that are related 
to the formation of MSPs: encounter rate and metallicity. This places an intimate relationship between 
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Figure 2: Gamma-ray and radio pulse profiles of PSR J1823-3021A in the globular cluster NGC 6624 jFreire et allboilh . 
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the observed 7-ray emission from globular clusters and the MSP population. The new findin gs that optical 
and i nfrared photon energy densities are also correlated to L 1 conform with the IC model (jCheng et al 



20101 ) . in which both the number of MSPs in a globular clu ster and the soft photon energy densities play 

(2011) to combine two of the above 



Hui et al 



important role in the resulted 7-ray output. This prompted[ 
cluster parameters and to carry out a 2-dimensional regression analysis. The best-fit relations are given 
by 

logL 7 = (34.12 ±0.29) + (0.42 ± 0.17) lo g r c 
logi 7 = (34.70 ±0.30) + (0.39 ±0.18) logT c 
logi 7 = (35.21 ± 0.29) + (0.49 ± 0.18) [Fe/H] 
logi 7 = (35.61 ±0.16) + (0, 



(0.62 ±0.29) log u opt 
- (0.96 ±0.49) log u m 
-(0.44 ±0.31) log u opt 



l ±0.17) [Fe/H] + (0.76 ±0.50) log u m 
The edge-on views of the best-fit fundamental-plane relationships are depicted in Fig. [31 



(1) 
(2) 
(3) 
(4) 



As discussed, the encounter rate is related to the number of dyn amically- forme d binaries while metal- 
licity may represent the likelihood of forming MSP in such binary (|Ivanoval 120061) . Therefore we expect 
that the number of MSPs in a globular cluster increases as the encounter rate and metallicity in the clus- 
ter. Furthermore, soft photons act as seed photons to be up-scattered to 7-ray energies by the relativistic 
particles originated in MSPs or their surroundings in the IC models. We here attempt to perform a 
3-dimensional regression analysis using encounter rate, metallicity, and one of the two soft photon energy 
densities: 

log L 1 — A log r c + B [Fc/H] + C log u m + D (5) 




Figure 4: Best-fit fundamental planes using a 3-dimensional regression analysis. Stellar encounter rate, metallicity, and one 
of the two soft photon energy densities are used. 



Table 2: Best-fit parameters in the regression analysis 



Equation© Equation© 



Parameter 


Best-fit value 


Parameter 


Best-fit value 


A 


0.22±0.15 


A 


0.23±0.14 


B 


0.35±0.14 


B 


0.34±0.13 


C 


0.51±0.41 


C 


0.41±0.24 


D 


35.14±0.31 


D 


34.78±0.34 



log L 7 = A log r c + B [Fe/H] + C log u opt + D (6) 

The best-fit results are tabulated in Table [2] and the edge-on views on these fundamental- plane rela- 
tionships are shown in Fig. 2J 



4. Models for 7-ray emission from globular clusters 



Two kinds of models for 7-ray emission from globular clusters have been widely discussed in the 
literature. In pulsar magnetosphere models, 7- ray emission is produced via curvature radiati on within the 



2008 



Venter et al 



2009). The pulsed 



magnetospheres of MSPs in the clusters (e.g.. IVenter fe de Jager , 
emission of a number of MSPs is superposed at different frequencies such that unless a small number 
of luminous MSPs dominate the combined 7-ray emission, as is the case for NGC 6624, pulsations from 
individual pulsars cannot be detected. 

In another class of models that involve inverse Compton (IC) emission, electrons are accelerated 
close to the MSPs or (re-)accelerated i n colliding wind shocks and scatter off the o ptical, infrared, or 



2007 



Cheng et all |201fj) , giving rise to the 



cosmic microwave background photons (jBednarek fe Sitareld . 
observed GeV emission. In this scenario, the 7-rays are intrinsically unpulsed. 

The 7-ray spectrum and especially the cutoff energy as observed for a number of GCs are generally 
similar to nearby MSPs in the Galactic field. Such coincidence may suggest that the observed 7-rays are 
curvature radiation from pulsar magnetospheres. On the other hand, averaged spectra may not be 
good discrimi nators in test i ng va rious models, because the I C models can al s o expl ain the GeV spectra 



equally well (jCheng et all l201Clh . Moreover, according to ICheng fc Taaml (|20031) . MSPs in globular 



clusters may be very different from those in the Galactic field based on radio and X-ray properties of 
MSPs in globular clusters. In fact, the former may^possess complicated magnetic fields due to frequent 



stellar encounter (jCheng fc Taaml . 120031) . thereby strongly affecting the polar/outer gap structures and 
in particular quench the outer gap. This s cenario is supported by the fact that a majority of the MSPs 
in 47 Tucanae are thermal X-ray emitters ([Bogdanov et al.1 . 120061 ) . 

The differences in the emission signatures of the two models should provide a diagnostic. In the pulsar 
magnetosphere model, 7-rays up to a few GeV is emitted, while IC processes may give rise t o 7-rays up 



to Te V energies. A candidate globular cluster reported to emit >10 GeV emission is Liller 1 ([Tarn et al 



20111 ). LAT observations above ~10 GeV as well as VHE observations by Cherenkov telescopes are 



crucial in testing both classes of models. 

The size of the emission region also differs between the two classes of models. The 7-ray emission from 
pulsar magnetospheres is more compact and that from the IC model is more extended. The observed 
7-ray emission region is consistent with a point source for almost all globular clusters and shows no sign 



of extended emission based on 1.5 year of data (Abdoetal 



2010bTI . While it is ex pected in the pulsar 

2010l) 



Cheng et al 



magnetosphere models, IC models predict a certain degree of extension. For example 
suggests that the IC emission size should be >10 pc, which in the case of 4 7 Tuca nae corresponds to >8' 

(|2010bh . Therefo re, unless the 



Abdo et al 



Cheng et al 



(see also their Fig. 4), larger than the upper limit imposed in 

actual diffusion coefficients are lower than that assumed in their calculations, the model of 
f|2010h may not be able to explain the compactness of the emission size. 

Current angular resolution of LAT is insufficient to resolve any extended emission significantly smaller 
than 10'. Observations at energies >10 GeV may better constrain the emission size in the future. This 
may involve continuous LAT survey observations and targeted observations by current and planned 
Ch erenkov telescopes, su ch as H.E.S.S.-II, MAGIC-II, VERITAS, and CTA. 



Clapson et al.l ([20111 ) present radio data taken using Effelsberg 100-m telescope in the vicinity of 



Terzan 5 which may impose some constraints on the IC models. The measured radio flux at 11 cm in 
the circle of radius 0.15° around the core of Terzan 5 (region 1 in their Fig. 1) is (1.41±0. 21) Jy (local 



backg round emission contributes about 30-40% of this flux). According to Eq. (27) of ICheng et al 



(|2010l ). the modeled radio flux at 11 cm is about 10 Jy and 6.8 Jy within a few arc-minutes when the 
seed photons in their IC model is infrared and optical photons, respectively. It is also unlikely that 
the magnetic field is much lower than 10 -6 Gauss used in their Eq. (27). It may still be possible that 
low-energy electrons (that is responsible for sy nchrotron radiation below the peak at ^44 GHz) diffuse 
further out to region 11 in lClapson et al.l (|201ll) 's Fig. 1, thereby explaining the enhanced radio emission 
in region 11, whose flux is (3.86±0.34) Jy, while high-energy electrons having a much shorter diffusion 
length stay close to the core upscatter ambient photons and give rise to the GeV emission up to 10 GeV. 
Detailed modeling is needed to understand the astrophysical conditions there. 

Recentl y, 7-ray emission from particles accelerated by non-accreting white dwarfs (WDs) was also 
suggested (jBednarekL 120121 ). Due to the large number of WDs, these particles may also contribute 
to the observed 7-rays from globular clusters. In this scenario, particles are accelerated in the inner 
magnetosphere of WDs, diffuse out, and produce 7-rays via IC scattering off various soft photon fields, 
e.g. star light, infrared radiation, and cosmic microwave background radiation (CMBR). 



5. TeV observations of globular clusters 



Gamma-ray emission above 100 GeV can be produced via Inverse Compton s cattering off variou s 
soft photon fields by relativistic leptons accelerated in the pulsar magnetosphere (jVenter et all 120091 ) . 
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or leptons re- accelerated in shocks generated by pulsar winds (jBednarek fc Sitarekl 120071) . The spec- 
trum of >100 GeV emission differs between various models, as illustrated in Fig. [SJ These expectations 
prompted the observations of a couple of globular clusters (i.e., 47 Tucanae, M5, M13, M15, NGC 6388, 



above 100 GeV but they have not resulted in any detection ( Kabuki et al. , 


2007; 


Aharonian et al.. 


2009; 


Abramowski et al., 


2011a: 


Anderhub et al.. 


2009; 


McCutcheon et al. 


. 2009). 



Recently , the H.E.S.S. collaboration announced the 7-er detection of >0.4 TeV emission in the vicinity 
of Terzan 5 (jAbramowski et all 1201 lbh . They estimated the probability of chance coincidence of this new 
VHE source, HESS J1747-248, with an unknown AGN or PWN to be 10~ 4 . The reported photon flux 
of 1.2xl0 -12 cm _2 s _1 between 0.44-24 TeV corresponds to 1.5% of the flux of the Crab nebula in this 
energy range. The spectrum can be fit by a simple power law of photon index T 7 = 2.5 ± 0.3 s t a t ± 0.2 sys . 
Unlike the 0.5-20 GeV emission, the TeV emission region is offset from the cluster core (at 2-er level) 
and shows sign of extended feature. These morphological features are hard to be reconciled in models 
where emission is produced directly, either magnetospheric or IC, by the MSP population in Terzan 5. 
Interestingly, radio emission at 11 c m from the Effelsbe rg Galactic plane surveys was detected 0.8° 
towards the north-west of Terzan 5 (jClapson et all 120111) . The authors suggested a PWN scenario in 
which rclativistic leptons are accelerated by the pulsar population, diffuse out, and produce the radio 
emission and TeV emission by synchrotron and IC emission, respectively . This scenario is supported 
by the diffuse, non-thermal X-ray emission detected earlier by lEger et all (|2010l ). A large drawback in 
this picture is that it cannot explain the different offset directions of radio(north-west)/TeV(south-east) 
emission from the cluster core. A recent st udy did not reveal any diffuse X-rays from six 7-ray emitting 
globular clusters (|Eger fc Domainkol 120121 ) . 

Another leptonic scenario that we consider here is that the TeV emission is originated from IC emission 
by electrons accelerated in colliding shocks between collective pulsar wind from the cluster and the 
Galactic wind, which can explain the fact that the TeV emission site is on the way between the cluster 
core and the Galactic centre. 

In leptonic scenarios, as noted by lAbramowski et all (|2011bl ) and iDomainkol (|201lh . there should be 
a Klein-Nishina cutoff at a few TeV, which is not observed. However, the lack of a cut-off may simply be 
a r esult of low photon statisti cs at these high energies. 



Abramowski et al 



(|2011b| ) further discuss two hadronic scenarios: cosmic-rays accelerated by a past 
supernova or in a short GRB remnant. They have the advantage of being able to explain the observed 
simple power law spectrum up to 20 TeV. In the supernova interpretation, they argue that given the 
lack of molecular clouds and thus the low interstellar medium density at this location (n m 0.1 cm~ 3 
is assumed), the cosmic-ray energy needed to produce the observed TeV flux reaches 10 51 erg, which 
is rather high for a supernova. The more exotic short GRB scenario is further developed in 



Domainko 



(|201ll ) , but a similar energy is required and it is unclear why this particular short GRB can b e so energetic 



Nakar, 



2007) 



(i.e., at least 10 51 erg, which is at the high end of isotropic-equivalent energy of short GRBs; 
and efficient in transferring the energy to cosmic-ray particles. In this scenario, thermal X-rays at the 
flux level of 10 -12 erg cm _2 s _1 are produced in sho cks caused by sub-relativist ic ejecta expelled during 
the merger event that heat the interstellar medium (jDomainko fc Ruffertl . [2008) . 

In summary, the origin of HESS J1747— 248 is unclear. Future X-ray observations at the centre of the 
TeV emission, as well as deeper TeV observations, are crucial in tackling the origin of the TeV emission. 
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2 -1 
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Figu re 5: Model compa rison for Terzan 5. The green lines rep resent the curvatur e radiation spectrum (and errors) 
from IVenter et_al. The m agenta lines are IC spectra from ICheng et all fcnioh . The black line shows an IC pre- 

diction from [ lie dnarek fc Sitarekl j2007l ). The bl ue lines represent I C spectra for different values of magnetic field in the 



I hoiih . 



Data arc taken from 



observations. The figure is taken from I Venter et aL 1 2011 ) 



Abdo et all l l2010bh . The red line shows the H.E.S.S. sensitivity in a 50-h 
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6. Gravitational waves from globular clusters? 



Double neutron-s tar (NS) mergers are believed to be progenitors of short gamma-ray bursts (GRBs; 



Naravan et all 1 19921 ) . It is expected that double NS mergers are major sources of gravitational waves, 



especiallyshortly before the final coalescence. Gravitational wave detectors, including LIGC0, Virgc^l, 
GEO-HF0, and LCGT§|, are yet to detect the first s ignal from thes e events. In our Galaxy, there are about 
ten confirmed and candidate double NS binaries (jLorimeii 120081 ). By adding the estimated coalescence 
rates derived from each known double NS binary system in our Galaxy, their merger rate is estimated to 
be USt^M yr" 1 . In t his case, the detection rate of advanced LIGO of extragalactic mergers may reach 



265Z?78 Y r (|Lorimen . 120081) 



Based o n spatial offset distri bution of short GRB progenitors from their host galaxies. iGrindlav et al 



( 20061 ) and IChurch et al.l (|2011f) have suggested that at least a portion of short GRBs arise from bina- 
ries formed dynamically in globular clusters. Of the ten Galactic double NS binaries (or candidates), 
PSR B2127+11C, itself a dynamical l y form ed one, is the only one situated in a globular cluster, M 15, a 
7-ray source candidate (jAbdo et al.l . l2010bl ). Due to the high encounter rate and high concentration of 
binaries in globular clusters, there should be a much larger population of double NS binaries in globular 
clusters that have escaped timing measurements. Therefore, the actual merger rate of Galactic and ex- 
tragalactic double NS binaries may be even higher, thereby increasing the chance to detect gravitational 
waves from globular clusters. 
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